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Abstract 

We calculate the two-loop QCD corrections to paracharmonium decays r/ c — * 77 and rj c — * gg 
involving light-by- light scattering diagrams with light quark loops. Artificial large mass expansion 
and convergence improvement techniques are used to evaluate these corrections. The obtained 
corrections to the decays n c — ► 77 and r] c — ► gg account for —1.25% and —0.73% of the leading 
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I. INTRODUCTION 



Since the year 1974, when bound cc states (charmonium) were discovered, they have been 
thoroughly studied experimentally and theoretically, contributing to better knowledge of the 
Standard Model parameters. Generally, the theoretical description of charmonium decays 
requires the knowledge of the perturbative effects in cc annihilation, and non-perturbative 
corrections to the charmonium wave function. However, ratios of decay rates into different 
final states are only sensitive to perturbative QCD contributions. In the present paper we 
focus on the subset of two-loop perturbative corrections to the ground-state paracharmonium 
decays rj c — > 77 and r\ c — > gg. Abelian 0(a s ) corrections to both decays are related to the 
one-loop correction to parapositronium decay by a trivial replacement of coupling constants 
fl. Non-abehan one-loop contributions were presented in Q, 3- 0(a 2 s ) corrections to 
Vc 77 were discussed in but those results did not include light-by-light scattering type 
contributions. Then we calculate these diagrams which form a gauge- independent subset. 

Multiloop corrections in the perturbation theory can be calculated by the asymptotic 
expansion [5|. But at the present stage, it is impossible to calculate the two-loop diagrams 
considered in this paper by the ordinary asymptotic expansions because of the complexity of 
the diagrams. In order to calculate them we adopt the method of the large mass expansion 
(LME) Q, IT!. We introduce an artificial large mass M for the internal charm quarks and 
execute the asymptotic expansion with the mass scale of the external charm quarks m c (soft 
scale) and the mass scale of the internal charm quarks M (hard scale). When the loop 
momenta of the internal charm quarks are factorized into the soft scale, we can expand the 
propagators of the charm quarks by the ratio m c /M. The expansions of the propagators 
reduce the original diagrams to the simpler ones which we can calculate. After that we 
obtain the results in a series of m c /M and identify the artificial mass M with the original 
one m c . 

The present paper is organized as follow. In Sec. Ull we review the example of para- 
positronium decay p-Ps — > 77 to demonstrate the method of LME and obtain two-loop 
QED results. In Sec. 1111} we adapt the results in the p-Ps decay to the paracharmonium 
decays r\ c — > 77 and r\ c — > gg, and obtain the two-loop QCD corrections by counting the 
color factors. In Sec. IIVI we have a summary. 
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FIG. 1: Two diagrams which contribute to the p-Ps decay at tree level. 
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FIG. 2: Two diagrams M\ ree and Mfp e are shown. 



II. PARAPOSITRONIUM DECAY 



a. Tree-level decay rate Two diagrams in Fig. [T] contribute to the tree level decay. 
If we artificially set the mass of the internal propagator M ^> m (with electron mass m), 
the total e + e~-annihilation cross section and the p-Ps decay rate are calculated in the low 
momentum limit as 



a (e + e 



> 77) 

r(p-Ps — > 77) 
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2(3 (1 + r 2 ) 2 ' 

(2(3) -Aa (e + e-^ 77) • 1^(0) | 2 , 



(1) 
(2) 



where (3 is the electron velocity, r = at e /m, and r = m/M. Here the value of the positronium 
wave function at the origin is ^(O)] 2 = l/7r(2a ) 3 with a = l/(ma e ). Making r = 1, we 
obtain the well known p-Ps decay rate at the tree level. The same result can be obtained 
using the optical theorem: 
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where two diagrams M^ ree and M^ 66 are shown in Fig. [21 However, we here expand the 
Wick rotated propagator by a large mass M before the loop integral as l/(/ 2 + M 2 ) = 
Er=oH 2 )7 M2(n+1) - The ri § ht hand side of E( l- © is calculated as 



Im[Mf ee + M^ ree ] = vra 2 (8r 4 - 16r b + •■■), 



(4) 



which coincides with the first terms in the expansion of Eq. ([T]). We introduce the gauge pa- 
rameter in the photon propagator to check the gauge symmetry. In general the introduction 



iMi = 




FIG. 3: Three diagrams Mi,M2 and M3 are shown. 

of the artificial mass M breaks the gauge symmetry, namely, the Ward-Takahashi identity. 
However, when both of electron and positron rest, the gauge symmetry is restored and the 
result in Eq. PJ) does not depend on gauge parameter. 

b. Two-loop correction We proceed to two-loop contributions involving light-by-light 
scattering diagrams. We calculate the corrections Act to the leading order cross section a 
through the optical theorem: 
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ImfMi + Ms + Mg], 



(5) 



2(3m 2 

where the three- loop diagrams Mi, M2, and M3 are shown in Fig. [31 Here the internal 
fermion loop is considered massless, while the case of electron in the internal loop has been 
considered in [8[]. The sum of imaginary parts of these diagrams is gauge- independent. For 
divergences, we employ the dimensional regularization with D = 4 — 2e. The diagrams with 
inverse orientation of the fermion loop exist and their contributions are same. Assigning the 
artificial large mass to the two internal electron propagators, we now have to consider several 
possibilities for the loop momentum scales. In the case of diagram M 1; for example, we take 
Ps,P7, and P12 as the three independent loop momenta and the following four regions exist : 1. 
p 3 < M, p 7 < M, p l2 < M, 2. p 3 ~ M, p 7 < M, p 12 < M, 3. p 3 ~ M, p 7 ~ M, p 12 < M, 
and 4. p 3 ~ M, p 7 <C M, pi% ~ M. To calculate the hardest three-loop integrals in region 
1 (all soft momenta), we use the MINCER package [9( and the factorized integrals are 
calculated in FORM 1 101 using well-known formulae. Although each of Im[Mi] and Impvy 
has the divergences of 1/e and depends on gauge parameter, the sum is free from divergences 
and gauge parameter. We can obtain it up to the first eight terms: 



Im[Mi + M 2 ] 



7T 



a' 



+r 4 (-0.394 + 41ogr) +r 6 (0.121 - 81ogr) 
4 



+ r 8 (0.113 + 12.3 log r) + r 10 (-0.316 - 16.9 logr) 

+ r 12 (0.481 + 21.7 logr) + r 14 (-0.618 - 26.7 logr) 

+ r 16 (0.734 + 31.9 logr) + r 18 (-0.834 - 37.3 logr). (6) 

Here and hereafter only three digits in the coefficients are shown. Im[M 3 ] can be calculated 
up to the first six terms: 

Im[M 3 ] • — A = + r 4 (-0.313 + 81ogr)+r 6 (-0.706 - 161ogr) 

a 

+ r 8 (1.40 + 25.0 logr) + r 10 (-1.93 - 34.8 logr) 

+ r 12 (2.31 + 45.2 logr) + r 14 (-2.62 - 56.1 logr). (7) 

Im[M3] is gauge invariant in itself and has no divergence. The series in Eqs. ([6]) and ([7]) 
apparently do not converge at r = 1. However, we can choose an expansion parameter z 
instead of r such that these series converge. The comparison between Eqs. (pQ) and (jl]) 
suggests an appropriate substitution, 

2r 2 



1 + r 2 ' 

The change of the variable reconstructs the series in Eq. (J6j) into the following series of z 



(8) 



ImfMx + M 2 ] ■ — = +z 2 (-0.445 + 0.51og^) + ^ 3 (0.166) 
+2 4 (0.0108 + 0.0111 log z) + 2 5 (-0.000899 + 0.00740 log z) 
+^(-0.00156 + 0.00462 log z) + z 7 (-0.00110 + 0.00297 log z) 
+/(-0.000690 + 0.00198 log^) + /(-0.000421 + 0.001371og^), (9) 

which are good convergent series. The series in Eq. ([7]) are also converted to good convergent 
series of z. Since the obtained series are limited in the first several terms, we should estimate 
the errors with the remaining infinite series. We rewrite the series in Eq. (Q as ^°^ 2 CiZ% 



and estimate the errors as 



Si^io Ci < l c 9l Si^i t % ' Here we can conservatively set t = 0.62 



because of the ratio c 9 /c 8 — 0.610. We can obtain the result, 



HM 1 + M 2l .-^(gc, + _^) ± J^ (10) 
= -0.2726 ±0.00034 (= d 12 ), (11) 
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FIG. 4: Diagram corresponding to Mi in Fig. [3] for the decay r/ c — > 77 (r/ c — > gg) is shown in the 
left (right). 

where we write this result as d\2 for the latter references. In the same way, we can obtain 
the results, 



Im[M 3 ] • — A = -0.4692 ± 0.0029 (= d 3 ). 
or 

ImfMi + M 2 + M 3 ] ■ — A = -0.7419 ± 0.0033 (= d tot ). 
or 



(12) 
(13) 



We check that the results in Eqs. ffTT]) and f[T2l are supported by the results of two summa- 
tion methods, Holder summation and Shanks transformation. 

III. CHARMONIUM DECAY 



c. Paracharmonium decays to photons The decay rate is given in the analogous form 
to Eq. © as 



rfac - 77) = m ■ ^a tot (cc - 77) • l^ c (o)p 



(14) 



where ip c is the charmonium wave function. The tree- level decay rate is To^c - > 77) = 
47rag(Q c ) 4 A' r c |'0 c (O)| 2 /m^ with the charm mass m c , the electric charge of the charm quark 
Q c , and N c = 3. We here use the color singlet wave function for r] c , 
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\Vc) 



Sij |cV), 



(15) 



where i and j are the color indices. Two-loop QCD corrections have the three diagrams 
corresponding to Mi, M2 and M3 in Fig. [3j Diagram Mi (cc — > cc) corresponding to Mi is 
shown in the left in Fig. HJ where t a is the generators of QCD. We restrict the quarks running 
in the fermion loop to the up, down, and strange quarks which can be regarded as massless. 
We can obtain Im[Mi (cc — > cc)}, counting the extra factors relative to Im[M 1 (e~e + — > e~e + )], 
as 



Im[Mi (cc -> cc)] 



Im[Mi(e-e H 



e~e i 



(Q c ) 2 [(Q„) 2 + (Q,) 2 + (Q s ) 2 ]-T^(iV c 2 



(16) 



where a s is the strong coupling constant and Tr = 1/2. The extracted factors in Eq. ([16]) 
are common to the diagram where the orientation of the fermion loop is reversed in Fig. H] 
and they are also common to the other two diagrams corresponding to M 2 and M 3 in Fig. 
Then, using Eqs. (JSJ), (|14l) and ([TBI) with Q u = Q c = 2/3 and Qd = Q s = — 1/3, we obtain 
the decay rate including the two-loop corrections as 

T( Vc -> 77) = r ( Vc -> 77) [l + (^) d tot \ . (17) 

d. Paracharmonium decays to gluons The decay rate is obtained by the formula anal- 
ogous to Eq. (jllj) . The tree-level decay rate is r (?7 c — > gg) = 87rafT^Cir| , c (O) | 2 / m c with 
Cp = 4/3. Three diagrams which correspond to three in Fig. [3] contribute to the two-loop 
QCD corrections. Diagram Mi (cc — > cc) is the right one in Fig. HI Using the wave function 
in Eq. (TT5T) . we obtain the common extra factors for Mi i2 (e _ e + — > e~e + ) as 

Im[M lj2 (cc -> cc)] = Im[M li2 (e-e+ -> e~e + )] ■ T 2 R C 2 F n q , (18) 

where n g is the number of the light quarks. The color factor for M 3 is different from one for 
Mi and M 2 , as T 2 r Cf(Cf — CA/^)n q with Ca = 3. Thus, we obtain the decay rate as 

r(Vc - 99) = r ( Vc ^ gg) [l + (^f ^ (c F d tot - \c A d^j\ . (19) 

It should be noted that if we average the colors of the charm and anti-charm quarks in this 
calculation, the divergences are not cancelled because the color factors for Mi and M 2 are 
different as Tr[t a tH c t d ]Tr[t a t d t c t b ] ^ Tr [t a t b t c t d ]Tr[t a t d t b t c ]. The QCD bare lagrangian has the 
gluon four-point interaction, unlike QED, and the gluon four-point function has divergences 
proportional to the interaction term in the bare lagrangian. However, the requirement that 
T) c is color singlet state in Eq. (fTBT) reduces the color factors for Mi and M 2 from the different 
ones to the common ones in Eq. ( 1181) and eliminates the potential divergences in the two-loop 
QCD corrections. 

e. Numerical values of corrections We first estimate the corrections in r\ c — > 77 as 

r(r? c -> 77) = r ( Vc - 77) [1 + (^)* 7 W + (2°) 

= r (77 c ^7 7 )[l + (-0.439) + (-0.0125 + 0.017 ■ 5°J hers )\. (21) 
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Here we use the values of the one-loop corrections ~ —3.38 Q, [ll| and the two-loop 
corrections are written as 5^ = dtot + 5° ther from Eq. (fl7|) . S° ther denotes the other two- 
loop corrections. We also use the coupling constant at the scale of the charm quark mass, 



a s (m c ) = 0.41 [121 ] . We find that the obtained two-loop corrections account for —1.25% of 
the tree-level decay rate, neglecting the uncertainty of d to t in Eq. f|T3|) due to its smallness. 
We next estimate the corrections in t] c — > gg in the same way: 

nvc - gg) = r ( Vc -> <^)[i + (0.633) + (-0.00729 + 0.017 • 5f ers )}. (22) 

where in the analogue to Eq. (|20|) . the one- loop correction is 8jp ~ 4.87 The 
two-loop correction is 5 { g 2) = n q (C F d tot - C A d 3 /2)/2 + 6°*^ from Eq. @H). The obtained 
corrections account for —0.73% of the tree level. The ratio of the decay rates in Eq. (l2Til to 
one in Eq. (j22|) is free from the ambiguity of the wave function ip c (0). To compare the ratio 
to the experimental value, we need the other corrections S° thers and §° thers ^ which we leave 
for future work. 



IV. SUMMARY 

We calculate the two-loop QCD corrections to the decay modes i] c — > 77 and r\ c — > gg 
involving light-by-light scattering diagrams with light quark loops by the method of LME. 
The obtained series in LME are transformed into ones with good convergence properties 
by a change of the expansion parameter. The obtained corrections to the decay r] c — > 77 
account for —1.25% of the tree-level and the corrections to r\ c gg account for —0.73%. The 
artificially introduced large mass does not break the gauge symmetry in the kinematic case 
of both initial particles at rest because the results are independent of the gauge parameter. 
Regarding the divergences of the two-loop QCD corrections to r] c — > gg in the present paper, 
the nature that the observed hadron, t] c in this case, is color singlet state, eliminates the 
potential divergences. 
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